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ABSTRACT 

We have investigated the impact of dissipationless minor galaxy mergers on the angular momentum of the remnant. Our simulations 
cover a range of initial orbital characteristics and the system consists of a massive galaxy with a bulge and disk merging with a much 
less massive (one-tenth or one-twentieth) gasless companion which has a variety of morphologies (disk- or elliptical-like) and central 
baryonic mass concentrations. During the process of merging, the orbital angular momentum is redistributed into the internal angular 
momentum of the final system; the internal angular momentum of the primary galaxy can increase or decrease depending on the 
relative orientation of the orbital spin vectors (direct or retrograde), while the initially non-rotating dark matter halo always gains 
angular momentum. The specific angular momentum of the stellar component always decreases independent of the orbital parameters 
or morphology of the satellite, the decrease in the rotation velocity of the primary galaxy is accompanied by a change in the anisotropy 
of the orbits, and the ratio of rotation speed to velocity dispersion of the merger remnant is lower than the initial value, not only due to 
an increase in the dispersion but also to the slowing -down of the disk rotation. We briefly discuss several astrophysical implications 
of these results, suggesting that minor mergers do not cause a "random walk" process of the angular momentum of the stellar disk 
component of galaxies, but rather a steady decrease. Minor mergers may play a role in producing the large scatter observed in the 
Tully-Fisher relation for SO galaxies, as well as in the increase of the velocity dispersion and the decrease in v/cr at large radii as 
observed in SO galaxies. 
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Numerical simulations, as well as observations, show that 
the final product of a merger between two disk galaxies de- 
pends mostly on their mass ratio. Major mergers of spi- 
ral galaxies, i.e. of pairs with mass ratios ranging from 
1:1 to 3:1, have been kno wn for decades to form pressure- 
supported elliptical galaxies dBarnes & Hernquistl 1991 ; Barnes ! 



tions. For example, they can induce changes that move late-type 
spiral galaxies towards earlier Hubble types, and contribute to 
the thickenin g and heating of the stellar disk 
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supported elliptical galaxies (tsarnes & Hernqu 
ll992tlBekki & Shiovalll997tlNaab et al.lll999r.l 



Bendo & Barnes 

200 0tlSpringel200aiCretton et al.ll2001t iNaab & Burkertll2003t 
iBournaud et alj2005l) . with "boxy" or "disky" shaped isophotes. 
Intermediate mass ratio mergers , with mass ratios from 4.5:1 
to 10:1, produce hybrid systems which typically have spiral- 
' like morphologies but elliptical-l ike kinematics (I Jog & Chitrd 
■ l2002t IBournaud et al.ll2004l 120051) . while minor mergers (mass 
ratio < 10:1) are much less violent dynamical processes 
that do not destroy the disk of the more massive galaxy 
in the merger. However, understanding the impact of minor 
mergers on galaxies is of particular interes t, considering that 
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can induce a number o f morphological signatures in the fi 
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The importance of minor mergers in shaping the evolution 
of galaxies has been supported by detailed numerical simula- 



WalkeretalJ 19 96: Velazquez & Whit dll999t Font et al 
iBenson et al.ll2004tlKazantzidis et alj |2008). The disks of galax- 
ies are generally not fully destroyed by these events, and a re- 
maining kinematically cold and thin component containing ~15 
to 25% of the initial stellar disk mass can indeed survive, e mbed- 
ded in a hotter and thicker disk (IVillalobos & Helmill2008L e.g.,), 
whereas the presence of gas in the disk g alaxy can largely pr e- 
vent the destruction of its thin stellar disk dMoster et al. 2009). 



However, while the role that minor mergers may play in the 
mass assembly, kinematics, and morphologies of present-day 
galaxies has been studied extensively both observationally and 
numerically, much less is known about the way angular momen- 
tum is redistributed during these events. Cosmological simula- 
tions show that there is a discrepancy between the specific angu- 
lar momentum of dark matter halos, which have undergone pre- 
dominantly minor mergers since z ~ 3, and the observed bary- 
onic component of present-day galaxies. The dark matter halos 
in cosmological simulations have spin parameters which are a 
few times smal ler than observed in disks of late-type bulgeless 
dwarf galaxies (Ivan den Bosch et al hOOltlD'Onghia & BurkerJ 
|2004|) . The evolution of the specific angular momentum of dark 
matter halos as caused by major and minor mergers can be de- 
scribed as a random walk process, with major mergers contribut- 
ing to the spin-up of dark halos, and with min or interactions 
mostly slowing them down dVitvitska et al.ll2002l) . 
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The consequence of an interaction is to transfer orbital an- 
gular momentum into internal rotation: at the beginning of the 
interaction, the most extended components first interact tidally, 
while the more tightly bound components feel the effect s of the 
encounter only in the final phases of th e merging process (iBarnesI 
11993 15 arnes in Kennicutt et al 1 fl998h . However, little is known 
about the detailed redistribution of the angular momentum be- 
tween the baryonic and dark matter components. iMcMillan et all 
d2007h have shown that, during major mergers of disk galaxies, 
the dark halos acquire only a small amount of rotation, with the 
dark component of the final remnant always having a stellar ro- 
tation speed to vel ocity dispersion, or v/cr, ratio less than 0.2. 
Ijesseit et al.l (120091) have discussed the amount of specific angu- 
lar momentum found in the stellar component of major merger 
remnants, showing that 1 : 1 mergers cause the destruction of the 
ordered rotation of the progenitor disks, while a portion of the 
original rotation can be preserved during mergers with higher 
mass ratios. If a dissipative component (gas) is present in the 
progenitor disks, gravitational torques exerted on it by the stel- 
lar disk can extract a significant fraction of its initial angular 
moment um, causing strong gas inflo ws into the circum-nuclear 
regions dBarnes & Hernquistl Il996l) . Finally, angular momen- 
tum redistribution also occurs during major mergers between 
spheroidal dissipationless galaxies, leading to a transformation 
of the orbital angular momentum into the internal angular mo- 
mentum of both baryons and dark matter, which can be efficient 
enough to prod uce fast rotating (v/cr > 1) stellar halos, as re- 
cently shown bv lDi Matteo et al] d2009h . 

A complete and detailed picture of the angular momentum 
redistribution during minor mergers is still lacking, even though 
minor mergers are expected to be m uch more common than 
major mergers jFakh" ouri & Mai 120081) and therefore obviously 
play an important role in g alaxy mass assembly. For exa mple, 
iDebattista et al. I d2006l) and lSellwood & Debattista I d2006l) have 
suggested that minor mergers could transfer angular momentum 
to the dark halo through the action of a stellar bar in the disk. 

This paper aims to investigate the evolution of the angular 
momentum in minor mergers and its redistribution between the 
stellar and the dark matter components of the primary galaxy 
and its satellite. We will investigate this process through the use 
of dissipationless N-body simulations of minor (mass ratio 10:1 
and 20:1) mergers between a disk plus bulge primary galaxy 
and a satellite having either an elliptical or a disk morphology. 
Specifically, we will show that independent of the morphology 
of the satellite galaxy or the orbital parameters, minor mergers 
always cause a slowing down of the stellar disk, as well as a 
loss of specific angular momentum within the disk, and that this 
process contributes to moving galaxies towards earlier Hubbl e 
types and more slowly rotating systems dEmsellem et al.ll2007l) . 
The decrease of the specific angular momentum of the stellar 
component of the disk galaxy is accompanied by an increase in 
the specific angular momentum of both dark halos (of the pri- 
mary and of the satellite), as well as of the pressure-supported 
stellar component that was initially part of the satellite. 

2. Models and initial conditions 

We studied the coalescence of a massive SO galaxy and an 
elliptical having a mass ratio 10:1 or 20:1. The massive SO 
galaxy (hereafter called gSO) is composed of a spherical dark 
matter halo and a spherical stellar bulge, both of which are 
not rotating initially, and represented by Plummer spheres 
dBinnev & Tremainell 19871) with total masses Mh and Mb and 
core radii r# and rg, respectively. The stellar disk is represented 



by a Miyamoto-Nagai density profile dBinnev & Tremainel 
Il987l) with mass M*, and vertical and radial scale-lengths of 
h„ and a„, respectively (Table Qj). The adopted bulge to total 
baryonic mass ratio of our SO galaxy model, B/T of 1:5, is 
bas ed on the most recent e stimate for the ratios of SO galax- 
ies dLaurikainen et alj|2007j). which are up to a factor 3 lo wer 
than earlier estimates (e.g..lSimien & de Vaucouleursll 19861) . As 
the work of iLaurikainen et alj d2007l) is based on K e band (2.2 
//m) data, it represents more realistic values for the o l d stel- 
lar component than those by ISimien & de Vaucouleursl dl986l) 
based on B -band data (the dust extinction being more signifi- 
cant at these wavelengths). The fraction of dark to baryonic mat- 
ter inside R50, the radius containing half of the baryonic mat- 
ter, is 1 6%, which is in broad agreement with observational es- 
tima tes (|WilUams et al.l|2009|) . and the resulting rotation curve 
(see IChilin garian et al. 2009) shows a rapid rise in the cen- 
ter, followed by a declin e at larger radii, typical of many SOs 

dNoordermeer et al]|2007l) . 

This paper is part of a larger program dDi Matteo et al.l2009t) 
whose aim is to study systematically the impact of major and 
minor mergers on angular momentum redistribution in galaxies. 
We analyze here the simple case of the impact of minor merg- 
ers on early-type, gas poor, SO galaxies and we plan, in future 
work, to extend this study to the whole Hubble sequence, with a 
variety of gas fractions and bulge to disk ratios. Besides, starting 
with such simple simulations allows us to compare our results 
with a well defined class of objects in order to verify the appro- 
priateness of our simulations and to gain insight in possible evo- 
lutionary avenues for relatively simple dynamical systems like 
SOs. 

The satellite elliptical galaxy, which is either ten (dEO, dEOl, 
dEOh) or twenty (sEO) times less massive than the gSO galaxy, 
consists of spherical stellar and dark matter components, both 
modeled with Plummer profiles, and both not rotating initially. 
The density profile of the satellite galaxy is the same in all sim- 
ulations, and we only changed the central density of the bary- 
onic component to study how changing the density would affect 
the angular momentum redistribution during the encounter. We 
thus considered "reference" dwarf E0 (dEO) models, as well as 
galaxies having a central volume density 50% higher (dEOh) and 
lower (dEOl) than the reference dEO galaxy (see TableQ}. 

To study the dependence of the results on the morphology of 
the satellite galaxy, we ran twelve additional simulations, mod- 
eling mergers of the giant gSO galaxy with a disk satellite (here- 
after called dSO), whose total mass is ten times smaller. The 
dSO consists of a spherical dark matter halo, a spherical stellar 
bulge, both of which are not rotating, and a rotating stellar disk. 
The morphological parameters of the dSO galaxies are gi ven in 
Table Q] We refer the reader to IChilingarian et al.l d2009l) for a 
more extensive description of these initial models. 

All the initial models and their individual components were 
evolved in isolation for 1 Gyr, before starting the interaction. 
At this time, the disk has reached a stable configuration, as dis- 
cussed in the Appendix. The galaxies were initially placed at 
a relative distance of 100 kpc, with a variety of relative ve- 
locities, to simulate different orbits. The orbital initial condi- 
tions for the gian t -dwar f interactions are given in Table 9 of 
IChilingarian et alj d2009l) . and we report here only on the orbital 
parameters of the gSO-sEO simulations (Table|2]i. 

We chose a reference frame with its origin at the barycen- 
ter of the system and an x-y plane corresponding to the orbital 
plane. The spin vector of the gSO galaxy defines initially an an- 
gle i\ - 33° or ;'i = 60° w ith respect to the z-axis (see Fig. 3 in 
IChilingarian et alj d2009l) ). For all simulations involving a dSO 
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Table 1. Galaxy parameters for the initial models. The bulge and 
the halo are modeled as Plummer spheres, with characteristic 
masses Mb and Mh and characteristic radii rg and r#. M» is the 
mass of the stellar disk whose vertical and radial scale lengths 
are h* and a,, respectively. 





gso 


dEOl 


dEO 


dEOh 


dSO 


sEO 


M B [2.3 x 10 lJ M o ] 


10 
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3.5 


M H [2.3 x 10 9 M o ] 


50 
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1.5 


M, [2.3 x 10 9 M o ] 


40 








4 




r B [kpc] 


2 


1.7 


1.3 


1.1 


0.6 


0.9 


r H [kpc] 


10 


2.2 


2.2 


2.2 


3.2 


1.55 


a, [kpc] 


4 








1.3 




h, [kpc] 


0.5 








0.16 





galaxy, we have chosen an initial angle between the dSO spin 
and the z-axis of 12 = 130°. The orbital angular momentum can 
be parallel (direct orbit) or anti-parallel (retrograde orbit) to the 
z-axis of the reference frame. 

When referring to specific encounters between the spiral 
and a satellite, the nomenclature adopted is a six(or seven)- 
character string: the first three characters are the type of the 
massive galaxy, always "gSO", followed by a three(or four)- 
character string for dEO (average concentration, dwarf ellipti- 
cal satellite), dEOl (low concentration dwarf elliptical), dEOh 
(high concentration dwarf elliptical), dSO (dwarf spiral satellite) 
or sEO (for the two times less massive elliptical satellite). This 
is followed by a string representing the type of encounter (see 
first column in Table [2] for the specific types), followed by the 
suffix "dir" or "ret", for direct or retrograde orbits, respectively. 
The final two numbers, which are either "33" or "60", indicate 
the initial inclination of the gSO galaxy with respect to the or- 
bital plane. For example, the nomenclature "gS0dE0h01dir33" 
refers to the encounter between the gSO galaxy (whose initial 
inclination i\ = 33°) and the high-concentration elliptical satel- 
lite dEOh, moving on a direct orbit, whose initial orbital pa- 
rameters are those corre sponding to id = Oldir in Table 9 of 
IChilingarian et al.l (120091) . 

All simulations (5 in total) were ru n usin g the Tree- 
SPH code described in ISemelin & Combesl d2002l) . A total of 
Ntot =528, 000 particles was used for all simulations, distributed 
between the gSO galaxy and the satellite (see Table [3). We also 
tested the dependence of the results on the number of particles 
used in the simulation, running some additional simulations with 
a total of Ntot- 1,056,000 particles. A Plummer potential was 
used to soften gravity on small scales, with constant softening 
lengths of e = 200pc (or e = HOpc for the high resolution 
simulations) for all particles. The equations of motion are inte- 
grated using a leapfrog algorithm with a fixed time step of 0.5 
Myr. With these choices, the relative error in the conservation of 
the total energy is close to 10~ 6 per time step. Since the work 
presented here only investigates simulated galaxies without any 
gas, only the part of the code evaluating the gravitational forces 
acting on the systems has been used. 

3. Results and discussion 

3.1. Angular momentum evolution 

As a result of the action of tidal torques and dynamical friction, 
the satellite loses its orbital angular momentum (hereafter AM) 
and gradually spirals deeper into the gravitational field of the 
primary galaxy, where it finally dissolves completely (meaning 
it does not survive as an entity within the core of the primary 



Table 2. Orbital parameters for the gSO-sEO interactions 



orbit id 


Ani 


1 mi 


L 


t, 


spin 


Oldir 


100. 


1.48 


29.66 


0. 


up 


Olret 


100. 


1.48 


29.66 


0. 


down 


02dir 


100. 


1.52 


29.69 


0.05 


up 


02ret 


100. 


1.52 


29.69 


0.05 


down 


03dir 


100. 


1.55 


29.72 


0.1 


up 


03ret 


100. 


1.55 


29.72 


0.1 


down 



a Initial distance between the two galaxies, in kpc. 

b Absolute value of the initial relative velocity, in units of 100 kms' 1 . 

c L =| ri„i x Vi„i I, in units of I0 2 kms^ 1 kpc. 

6 E = v,„, 2 /2 - G(m { + m 2 )/r ini , with mi = 2.3 x 10"M o and m 2 = 

2.3 x 10 10 M o , in units of 10 4 fcmV 2 . 
e Orbital spin, if parallel (up) or antiparallel (down) to the z-axis 



Table 3. Particle numbers for disk galaxy and satellites 



gSO 


dE01,dE0,dE0h,dS0 
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Fig. 2. Evolution of the absolute value of the internal angular 
momentum in minor mergers with different baryonic mass con- 
centrations of the merging satellite. The mass of the satellite is 
the same, only the central volume density is 50% higher (dotted 
lines) or 50% lower (solid lines) than in the reference satellite. 
The internal angular momentum of gSO galaxy and satellite are 
shown in red and blue, respectively. Left panel: direct mergers. 
Right panel: retrograde mergers. 

galaxy). In this section, we discuss how the orbital AM is con- 
verted into internal AM and how it is redistributed among the 
different components (stars and dark matter) of the two galaxies. 

3.1.1. Redistribution of orbital into internal angular 
momentum 

Fig. Q] shows the evolution of the absolute value of the total, or- 
bital and internal angular momenta of the system, for some of 
our simulations. In all the cases, the total (internal + orbital) AM 
is conserved during the interaction, with an accuracy of < 10~ 6 
per time step. The orbital AM is constant until the first peri- 
center passage of the satellite. At this time, dynamical friction 
and tidal torques act on the system, converting part of the or- 
bital AM into internal AM. This continues to occur during each 
successive passage, until the two galaxies finally merge, and the 
orbital AM has been completely converted into internal rotation 
of the merged system. Looking at the repartition of the internal 
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Fig. 1. Top panels: Evolution of the absolute value of the total (solid line), orbital (dotted line) and internal (dashed line) angular 
momentum for some of the simulated 10:1 mergers, for both the stellar and dark matter components together. Bottom panels: 
Evolution of the absolute value of the internal angular momentum of the massive gSO galaxy (red dashed line) and of the dEO satellite 
(blue dashed line), for both the stellar and dark matter components together. The angular momentum is in units of 2.3 x 10 1 l M G kpc 
km s . 
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Fig. 3. Evolution of the absolute value of the internal angular 
momentum of the stellar and dark matter components of the gSO 
galaxy (solid and dotted red lines) and of its satellite (solid and 
dotted blue lines) for four 10: 1 mergers. The angular momentum 
is in units of 2.3 x 10 1 l M Q kpc km s _1 . 



AM between the two galaxies (Fig. Q~| bottom panels), we can 
see that the satellite galaxy always gains part of the orbital AM, 
while the gSO galaxy either gains or loses angular momentum, 
depending on the relative orientation of its spin and of the orbital 
AM (direct and retrograde orbits, respectively). Note that since 
in retrograde encounters the orbital AM is anti-parallel to the 
z-axis of the reference frame, their total AM is systematically 
lower than that of the corresponding direct encounter. 

Changing the baryonic mass concentration of the satellite 
(Fig. |2]i affects the redistribution of the internal AM between 
the two galaxies, in the sense that, for a given orbit, the denser 



the satellite, the smaller the amount of AM it absorbs and the 
larger the variation in the internal AM of the primary galaxy. 
This trend is due to the fact that a satellite with higher density 
decays more rapidly in the central regions of the primary galaxy, 
since it is less susceptible to tidal effects of the primary galaxy, 
preserves a higher percentage of bound mass, suffers a stronger 
gravitational drag. Since tidal torques are less effective on denser 
satellites (see §3.2 in Di Matteo et al. 2009), the orbital angular 
momentum is more efficiently redistributed to the dark halo of 
the primary galaxy, the component which mainly drives its or- 
bital decay through dynamical friction. 

3.1 .2. Internal angular momentum of the stellar and dark 
matter components 

Here we investigate how the AM is redistributed in each of the 
two galaxies, between the stellar and the dark matter compo- 
nents. 

In Fig. [3] we show how the distribution of the internal AM 
between the baryons and dark matter evolves during the merg- 
ing process. Near the time of the first barycentric passage, the 
two dark matter halos (primary and satellite) acquire part of the 
orbital AM, leading to rotation of a halo that initially was not 
rotating. In general, the amount of internal AM absorbed by the 
dark matter component of the primary galaxy is greater (by a 
factor 2-3) for pairs on direct orbits than for those on retrograde 
ones. For retrograde orbits, we note that the dark halo of the 
merger remnant, rotates in the direction opposite to that of the 
stellar disk, because it absorbed part of the orbital AM. 

Also the stellar component of the satellite absorbs orbital 
AM (Fig. |3). For retrograde orbits, this conversion of orbital 
into internal AM is a mechanism for creating a counter-rotating 
old stellar component in a merger remn ant (as first noted by 
iKormendvll 1 984t iBalcells & Ouin for minor mergers in- 

volving elliptical galaxies). Note that, contrary to the simulations 
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bv lBalcells & Ouin nl (1199(H) . our satellite does not survive as a 
recognizable entity to the end of the merging process. Its distri- 
bution can be described as "thick disk-like", encompassing the 
remaining disk and extending out to ~10 kpc from the center of 
the remnant. Moreover, less than 4% of the internal stellar AM is 
transferred to the remnant's bulge by the end of the merger pro- 
cess. Minor mergers do not add significant amounts of angular 
momentum to the bulges of SO galaxies. 

Although none of the initially pressure supported compo- 
nents were rotating at the beginning of the simulation, they later 
acquire part of the orbital AM and begin to rotate. The case of 
the halo component is simple: the more massive and extended, 
the more AM it absorbs. However, the evolution of the AM of 
the stellar disk of the primary galaxy is quite different. In all 
simulated orbits, independent of the orbital parameters and the 
morphology of the satellite galaxy (whether pressure- or rota- 
tionally supported), the disk of the primary galaxy always loses 
AM during the interaction. Although the decrease of the internal 
AM of the disk varies depending on orbital parameters and on 
the internal structure of the satellite (see Fig. [3}, on average the 
relative change of the internal AM is AL/L~ 10% for mergers 
with a mass ratio 10:1. We also find that satellites with higher 
central stellar densities lead to remnants with a lower internal 
AM in the stellar component that was formerly part of the disk 
of the primary. 

3.2. The slowing-down of rotating disks 

The decrease of the internal AM of the primary stellar disk dis- 
cussed in the previous section affects the whole of the disk, up 
to distances of about 5 to 10 times /?5fEl To better demonstrate 
the loss of angular momentum, in Fig. [4] we show the evolu- 
tion of the specific internal AM, l m , -< E,r, x v, >, of the stars 
that were initially part of the primary galaxy. We have evaluated 
li nt for five different radial regions of the galaxy: r < O.SRso, 
0.5R 50 < r < R 5a , R 5 o < r < 2R 5Q , 2R 50 < r < 5R 5a , and 
5/?5o < r < 10/?5o. In all cases considered here, we find that the 
stellar component slows down out to at least 5/?so, i.e. the radius 
containing ~ 80% of the total stellar mass (Pig. [5j, regardless of 
the initial orbital parameters, inclination angles, or satellite mor- 
phologies. The same behavior is also found in the merger with 
the disky satellite, before it is ultimately destroyed by the tidal 
field of the primary galaxy. 

The question poses itself what this decrease in rotation speed 
is due to? During the merging process, the variation of the 
half mass radius, R^q, is not more than 10%, and what varia- 
tion there is in the mass distribution is not sufficient to explain 
the observed decrease. To gain further insight, we evaluated the 

anisotropy parameter, y3 = 1 - , where o> and cr, are the 
velocity dis persions of the radial and t angential components, re- 
spectively dBinnev & Tremainelll987h . If the velocity distribu- 
tion is isotropic then f3 = 0, and if the velocity distribution is 
dominated by tangential motions then f3 < 0, and fj > if 
dominated by radial motions. As shown in Fig. [6] the merging 
process is accompanied by an evolution of the /3 parameter in 
the stellar component. The change in /? is particularly large out- 
side of /?5(). In these outer regions, the stellar orbits, which be- 
fore the interaction were dominated by tangential orbits, tend 
to become increasingly radially dominated as the merger ad- 
vances. While the outermost region, 2Rso < r < 5Rsq, of the 
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Fig. 7. The evolution of radial and tangential fractional velocity 
dispersion in different radial regions: r < Q.SR^o, 0.5/?5o < r < 
Rso, Rso<r< 2R 50 , 2R 50 <r< 5R 50 and 5R 50 < r < 10R 50 . 



merger remnant is still dominated by tangential motions, the mo- 
tion of stars in the region R$o < r < 2R^o can become nearly 
isotropic due to the interaction with the satellite (see for exam- 
ple the gS0dE001dir33 simulation shown in Fig|6]i. A similar 
change in the ani sotropy was previousl y noted in simulations 
of s ingle mergers ( Bournaud et al. 2005b and o f multiple merg- 
ers (lAthanassoulal 2005 : Bournaud et al.ll2007l) . In general, the 
strongest evolution of the /3 parameter is found in direct encoun- 
ters. 

If the velocity dispersion of a disk increases in a merger, then 
the effective rotational velocity decreases, due to the asymmetric 
drift, see, e.g., §4 of Binney & Tremaine (1987). The anisotropy 
parameter is a way to represent the conversion of rotational mo- 
tion into random motions. To confirm that asymmetric drift is the 
underlying cause, we evaluated the evolution of the radial and 
tangential components of the velocity dispersions in our simu- 
lations. In FigfT] their fractional evolution (i.e., cr r (t)/o- r (t = 0) 
and cr,(f)/cr,(f = 0)) is shown for one of our simulations. 

The minor merger is accompanied by an increase of the ra- 
dial velocity dispersion at all radii, but particularly in the outer- 
most regions. At the same time, the transverse velocity disper- 
sion decreases throughout the whole disk, except in the inner re- 
gion (inside 0.5/?5o), where the contribution of bulge stars, which 
acquire AM during the interaction, leads to an increase of <r t . 



Rso = 3.3kpc is the initial half-mass radius of the baryonic compo- 



nent. 



3.2.1 . The effect of increasing the merger mass ratio to 20:1 

In the previous sections we have discussed the impact of 10:1 
mass ratio mergers on the redistribution of AM, focusing partic- 
ularly onto the cause of the decrease in the rotation speed of the 
primary stellar disk. In this section, we want to focus on the ef- 
fect that increasing the mass ratio in mergers, specifically from 
10:1 to 20:1, has on the evolution of the anisotropy parameter 
and the decrease of the specific AM of the disk (Fig. [8}. Also 
for 20:1 mass ratio mergers a decrease of the /3 parameter is 
observed (again especially in the outer regions), indicating that 
the velocity dispersion is increasing in the radial direction. The 
overall impact on the disk rotation is, of course, less pronounced 
than for the 10:1 merger simulations. In this case, we find that 
the average decrease in the internal AM of the stellar compo- 
nent is AL/L k 8%, which is about 30 - 40% lower than the 
value found for 10:1 mergers. However, this is an average value, 
whereas the actual decrease in any particular merger depends 
on both the morphology and the central density of the satellite 
galaxy. 



6 



Qu et al.: Slowing down of galaxy disks in dissipationless minor mergers 



60 
50 
40 

J 30 
20 
10 

60 
50 
40 

J 30 
20 
10 



LgS0dE001dir33 



JLgS0dE004dir33 



LgS0dE001 ret33 



JLgSOdE001dir60 



JLgS0dE004ret33 



<R 5( /2 
R 50 /2<i<R 50 
R 50 <r<2R 50 
2R 50 <r<5R 50 
5R»„<i < 10R,, 



JLgS0dS001dir33 



0.5 



1.0 1.5 2.0 

t(Gyr) 



2.5 



0.5 



1.0 1.5 2.0 
t(Gyr) 



2.5 



0.5 



1.0 1.5 2.0 
t(Gyr) 



2.5 
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Fig. 6. Evolution of the anisotropy parameter, /3, of the stars in the gSO galaxy for the six merger simulations whose specific internal 
AM evolutions have been shown in Fig. [4] direct an d retrograde merger simulations with a dEO or a dSO satellite. The anisotropy 
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3.3. Rotation speeds and v / cr 



The decrease in the specific angular momentum of the stellar 
component during a minor merger is obviously reflected in other 
dynamical properties of the merger remnant. As shown in Fig. [9] 
the minor merger reduces the rotation speed of the stellar com- 
ponent of the disk galaxy, and increases its velocity dispersion. 
To evaluate these quantities, we considered all stellar particles 1 
kpc above and below the meridional plane of the primary galaxy, 
including those that were initially part of the satellite galaxy. In 
the case of a 10:1 direct merger, the rotation speed at r — 2Rso 
decreases by about 25%, from the initial value of v = 200km/ s to 
about v = 150A77i/s. At the same time the velocity dispersion in- 
creases over the whole disk, e.g. by about 30% at r = 2R50. This 
leads to an overall decrease of the v/cr ratio over the entire extent 
of the disk, e.g. from 2.3 to 1.4 at r — 2R$o (Fig. [9). The decrease 
in the rotation speed and v/cr depends on the merger mass ratio: 
at the end of the 20: 1 merger the remnant shows higher values 
of both v and v/cr than the 10:1 merger. Finally, the decrease 
in the rotation speed in the remnants of the two retrograde sim- 
ulations shown in Fig. gS0dE001ret33 and gS0sE001ret33, 
is due to both a slowing down of stars in the gSO disk (as dis- 
cussed previously) and to a negative contribution coming from 
stars formerly in the satellite galaxy. These stars have acquired 
part of the orbital AM, which is in the direction opposite to the 
spin of the gSO, and they form a counter-rotating extended stel- 
lar component which contributes to the overall decrease in the 
line-of-sight velocity of the remnant galaxy, and to the increase 
in its velocity dispersion. 
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Fig. 8. Upper panels: Evolution of the absolute value of the spe- 
cific angular momentum of stars in the gSO galaxy undergoing a 
20: 1 direct (left panel) and retrograde (right panel) merger. The 
specific angular momentum has been measured in five different 
radial regions, just as in Fig. [4] Lower panels: Evolution of the 
anisotropy parameter,/?, for the mergers shown in the upper pan- 
els, and measured in the same five radial regions. 



4. Summary and Discussion 

In this paper, we have studied the redistribution of the angu- 
lar momentum during dissipationless minor mergers between a 
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Fig. 9. Upper panels: Line-of-sight velocities and velocity dis- 
persions of stars in the initial gSO model galaxy (dotted lines) 
and in the final remnant of a 10:1 (thick solid lines) and a 20:1 
(thin solid lines) merger. Direct encounters are shown in the left 
panel, retrograde encounters on the right. Lower panels: v/cr ra- 
tio of the initial gSO model (dotted lines) and of the final remnant 
of the 10:1 (thick solid lines) and 20:1 (thin solid lines) mergers 
shown in the upper panels. Distances are in multiples of R^. 

massive disk galaxy and an elliptical- or SO-like satellite with 
masses that are either ten or twenty times smaller than that of 
the primary galaxy. From an analysis of these simulations, we 
find that: 

- During the merging process the orbital angular momentum 
is redistributed into internal angular momentum of both the 
primary galaxy and the satellite; 

- The total internal angular momentum of the primary galaxy 
can either increase or decrease, depending on the orientation 
of the orbital spin (direct or retrograde orbit); 

- While the initially non-rotating dark matter halo of the pri- 
mary galaxy always acquires angular momentum, the spe- 
cific angular momentum of the stellar component always de- 
creases, independent of the initial orbital parameters or the 
morphology of the satellite galaxy; 

- The decrease in the rotation speed of the disk of the main 
galaxy is accompanied by a change in the distribution of 
the types of stellar orbits, especially outside of the half mass 
radius, R50. Generally, the radial component of the velocity 
dispersion becomes more and more important, as the merger 
progresses, thus increasing the anisotropy parameter,/?, from 
its initially negative value; 

- The ratio of the rotation speed to the velocity dispersion of 
the disk, v/cr, decreases at all radii due to both an increase 
in the velocity dispersion of the disk, which is heated during 
the merging process, and a decrease in the rotation speed of 
the disk. 

There are a number of astrophysical implications to the re- 
sults we presented of our (rather idealized) simulations, which 
we will now discuss very briefly. The Tully-Fisher, or T-F, re- 
lation between the rotation speed and total luminosity of disk 
galaxies has, among others, the following two interpretations. 
Either the T-F relation originates from the cosmological equiva- 
lence between t he dark halo mass and the circular velocity (e.g., 
iMo et al.ll 1998h . or it is due to self-regulated star formation in 



disks of different mass (e.g.. lSilkll997l) . The study of SO galaxies 
may help to constrain which of these hypotheses is correct, and 
thus to unravel the evolutionary history of SO galaxies. Studies 
of the differences between the Tully-Fisher relation of SOs and 
spirals, while showing some small discrepancies in the offset 
and scatter, generally find that the SO galaxies follow the spiral 
relationship, albeit w ith a la rger scatter (e.g.. iHinz et"aDl2003b 
lBedregaletai1l2006l but see IMathieu et all d2002l) ). The offset 
and larger scatter have been interpreted as revealing a generally 
older stellar population in SO galaxies than in spirals and a range 
of tim es during which the SO galaxies started to evolve pas- 
sively (IMathieu et al.l20"02h . However, as noted in lBedregal et al.l 
(120061) . there is no strong correlation between the age of the SO 
galaxies and the offset from the T-F relation for spiral galaxies. 
This could be due to an insufficient sample size, or it could indi- 
cate that other processes may play a role in the larger scatter in 
the T-F relation among SO galaxies. 

We have shown that minor mergers always lead to a decrease 
in the rotation speed of the primary galaxy. Decreasing its ro- 
tation speed while increasing its luminosity, will move an SO 
galaxy back towards the spiral T-F relation. While it is prema- 
ture to attempt to quantify this effect in the offset, it is clear that 
minor mergers would increase the scatter by counteracting, to 
some extent, passive evolution. In addition to slowing the rota- 
tion speed of the disk, we also found that the velocity dispersion 
of the disk increases, which decreases the v/cr of the disk, es- 
pecially in the outer parts. As we have show in Fig. [9] this may 
include an actual decrease with radius in the very outer parts 
of the disk. This is qualitatively similar to what is observed in 
the dynamic s of planetary nebulae in the outer disks/halos of 
SO galaxies dNoordermeer et aLll2008l) and it may suggest that, 
through simulating multiple minor m ergers, we may be ab le to 
reproduce the veloc ity structure (e. g .jBournaud et al1l2007l) . 

The models of IVitvitska et al.l (12002 ) have shown that, in 
mergers of dark matter halos, the specific angular momenta can 
either increase or decrease. Over time, this would lead to a in- 
crease in the dispersion of the specific angular momentum of the 
halos. Our results on the dark matter halo are consistent with this 
but we find systematically a decrease in the specific angular mo- 
mentum of the stellar component. So instead of a "random walk" 
over time with only an increase of the scatter of the angular mo- 
mentum, our results appear to suggest a systematic decrease with 
time in the angular momentum of the stellar component if minor 
mergers play an important role in the growth of galaxies and the 
redistribution of angular momentum. 

Finally, our results indicate that (single) minor mergers are 
capable of moving disk galaxies towards earlier types, since 
they reduce the specific angular momentum an d rotation speed, 
and in crease the velocity dispersion of the disk. iBournaud et al.l 
(I2007h have shown that a sequence of multiple minor mergers 
can lead to remnants having global morphological (flattening, 
Sersic index, etc.) and kinematical (v/cr ratio) similar to those 
observed in real elliptical galaxies. Within this context, it will be 
interesting to study if the angular momentum content of the stel- 
lar and dark matter components becomes more dispersion domi- 
nated and if size of this effect depends on the total merged mass 
rather than on the mass ratio of each merger, and if successive 
multiple mergers always slow down the stellar component of the 
primary galaxy. Many follow-up questions remain - for exam- 
ple, is the specific angular momentum of the stellar component, 
independent of the way mass has been assembled in major and 
multiple-minor mergers. What role does the dissipative compo- 
nent - gas - play in the way angular momentum is redistributed 
during the merger and in the remnant? How does the gas con- 
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tent affect the magnitude of the decrease in the rotation speed of 
the stellar disk? Does each minor merger contribute similarly to 
the net slowing down, or does each successive merger become 
less effective in changing the angular momentum and rotation 
speed? Does the net effect saturate, so that we will ultimately be 
left with nothing but very slow rotators? 

Our simulations show that a single dissipationless minor 
merger is sufficient to considerably reduce the v/cr ratio, typi- 
cally by a factor 1.5-2 at 2/?5o. Our merger remnants can still be 
classified as early-type disk galaxies, with a central bulge com- 
ponent and a disk that is thicker and kinematically hotter than 
that of the progenitor. If the cumulative effect of dissipationless 
multiple minor mergers leads to systematically slowing down 
stellar disks, then the final remnants, which should have an el- 
liptical morphology at that point, should not exhibit fast rotation 
at l arge radii. The existence of elliptica ls with fast rotating ha- 
los dRix et alJI 19991: ICoccato et al.ll2009h may thus require other 
formation mechanisms, such as maj or mergers of two ellipt ical 
progenitors, as recently suggested bv lDi Matteo et all (l2009h . or 
dissipationless mergers of dis k galaxies with higher mass ra- 
tios, ty_gicany_3d accord ing to lBendo & Barnesl ( 120001) ; but see 
also lCretton et alJ (l200lh . However, the presence of a dissipative 
component, whether in the progenitor disks or accreted after the 
merger, may change this picture. 

We will investigate all these points and their astrophysical 
implications as discussed above in subsequent papers. 
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Appendix A: Kinematics of the gSO galaxy when 
evolved in isolation 

The goal of this paper is to investigate the role minor mergers 
play in the angular momentum redistribution and kinematics of 
the massive galaxy. Therefore we must distinguish the impact 
of minor mergers on the properties of early type disk galaxies 
from the effects of secular evolution. To make this distinction, it 
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Fig. A.l. Upper panel: Line-of-sight velocities and velocity dis- 
persions of stars in the gSO galaxy which evolved in isolation, 
at t=0 (dotted lines) and at t=3 Gyr (solid lines). Lower panel: 
Corresponding v/cr ratios at these two times. Distances are in 
multiples of R50, 

Table A.l. v/cr ratios at different radii, for the gSO galaxy which 
evolved in isolation and after a minor 20:1 and 10:1 merger. The 
initial v/cr values are provided in Column 2. 







t=0 




t=3 Gyr 










isolated 


20:1 merger 


10:1 merger 


v/cr at r 


- ^50 


1.15 


1.03 


0.8 


0.72 


v/cr at r 


= 2^50 


2.44 


2.14 


1.57 


1.37 


v/cr at r 


= 4fl 50 


4.85 


3.73 


2.80 


1.81 



is also important to study the evolution of isolated galaxies. To 
this end, we have run a control simulation in which the massive 
SO galaxy is evolving as an isolated system over 3 Gyr, i.e., the 
same duration as the merger runs. 

As can be seen in Fig. Al, the kinematics of the isolated SO 
galaxy does not change significantly over this period: no increase 
is found in the velocity dispersion of the central regions or the 
outer disk, and also the final line-of-sight velocities are remark- 
ably similar to the initial values, apart from some insignificant 
decrease outside 2R50. Of particular interest is that the final v/cr 
ratio of the isolated massive galaxy differs significantly from that 
of the massive galaxy even after a minor merger, even in the 20: 1 
case. Comparing Fig.Al with Fig|9]shows that the minor interac- 
tion induces a decrease at all radii, even inside 2R$o, which is not 
found for the galaxy that evolved in isolation. Secular evolution 
causes only a slight decrease in v/cr - the final ratio is 10% lower 
than the initial one at r=/?5o, and 23% at r=47?go, whereas even a 
20: 1 merger induces a much more pronounced slowing down of 
the disk at all radii, with a final decrease in v/cr of 30% at r=Rso 
and 42% at r=4/?5o (see Table IA.ll . However, it will be of some 
interest to find out at what merger mass ratio the decrease in v/cr 
is comparable to that of secular evolution. 



